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00 ' We report the results of ab initio molecular dynamics simulations of liquid 

o : 

. Ga^^Asi-j; alloys at five different concentrations, at a temperature of 1600 K, 

. just above the melting point of GaAs. The liquid is predicted to be metallic 

I . at all concentrations between x = 0.2 and x = 0.8, with a weak resistivity 

-a ; 

^ , maximum near x = 0.5, consistent with the Faber-Ziman expression. The 

CJ ■ 

• • . electronic density of states is finite at the Fermi energy for all concentrations; 

^ . there is, however, a significant pseudogap especially in the As-rich samples. 

The Ga-rich density of states more closely resembles that of a free-electron 
metal. The partial structure factors show only a weak indication of chemical 
short-range order. There is also some residue of the covalent bonding found 
in the solid, which shows up in the bond-angle distribution functions of the 
liquid state. Finally, the atomic diffusion coefficients at 1600K are calculated 
to be 2.1 X 10"^ cm^/sec for Ga ions in Gao.8Aso.2 and 1.7 x 10~^ cm^/sec 
for As ions in Gao.2Aso.8- 
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I. INTRODUCTION 



The study of liquid metals and alloys has drawn considerable attention recently, in par- 
ticular due to the possibility of carrying out first principles or ab initio calculations for these 
systems. In these calculations the electronic structure is evaluated quantum mechanically 
using Density Functional Theory (DFT) and the corresponding forces are used to move the 
ions according to classical molecular dynamics. Using this approach it is possible to calcu- 
late both the atomic and electronic structure consistently and to see how changes in one are 
correlated with changes in the other. By now, a number of groups have used these methods 
to calculate both the thermodynamic and transport properties of a variety of liquid metals 
and alloys [HH^. 

Liquid semiconductors are of particular interest from the point of view of ab initio calcu- 
lations. By "liquid semiconductors," we mean liquids of materials which are semiconducting 
in their solid phases, such as Si, Ge, GaAs, and CdTe. Somewhat surprisingly, most of these 
are reasonably good metals in their liquid phases. For example. Si, Ge, and GaAs all have 
conductivities near melting which lie in the metallic range, and which tend to decrease with 
increasing temperatures, as is characteristic of metals. This metallic behavior is correlated 
with an increase in coordination number on melting, the liquid is thus more close-packed 
than the solid and has a higher density. By contrast, i-CdTe is poorly conducting in its 
liquid state and its conductivity increases with increasing temperature, characteristic of 
semiconductors. 

Recent ab initio calculations for several of these materials give behavior which is in good 
agreement with these experiments. Godlevsky et al have found, in agreement with ex- 
periment, that stoichiometric GaAs is metallic, whereas stoichiometric CdTe is a reasonable 
insulator. These differences in the electronic properties were related to the differences in 
the structural properties occurring within the melt. An earlier calculation by Zhang et al 
studied stoichiometric £-GaAs using the Car-Parrinello version of ab initio molecular dy- 
namics 0,0]. In this calculation too, it was found that ^-GaAs is a metallic, weakly ionic 
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liquid, with a larger coordination number than the insulating solid phase. 

In this paper, we describe a numerical study of £-Ga^Asi_j; over a range of concentrations, 
using ab initio techniques. The liquid is miscible over the whole concentration range, unlike 
the solid, which exists only at stoichiometry. Such a study is of interest for a variety of 
reasons. The properties of of i-Ga and ^-As stand in marked contrast to each other. l-Ga 
is a close-packed liquid metal with a coordination number of ~ 9 and its electronic density 
of states is almost free-electron like. In contrast, i-As has the same coordination number as 
in the crystalline phase (~ 3), and is a narrow band-gap semiconductor in the liquid state. 
Thus while we expect i-Ga. to show metallic bonding, i-As is expected to retain the covalent 
character of the bond upon melting. Previous ab initio calculations for these liquids 
have indeed confirmed this picture. Thus as the stoichiometry is varied for ^-Ga^-Asi-a;, we 
would expect interesting changes both structurally and electronically. It is also of interest 
to compare the structures of ^-Ga^jAsi-a; for x = 0.2 and x = 0.8 with the structure for 
the corresponding pure liquids (As and Ga respectively) to see how changes in structural 
properties are correlated with changes in electronic properties. 

We turn now to the body of the paper. In section II, we briefly summarize our approach 
and method of calculation. Our results are presented extensively in Section III, together 
with some analysis connecting the calculations to a qualitative picture of the liquid state. 
Finally, in Section IV, we give a short concluding discussion. 



II. METHOD AND COMPUTATIONAL DETAILS 

Our method of carrying out the ab initio simulations is similar to that of our previ- 
ous work P,p!0|. We use the plane- wave pseudopotential approach with generalized norm- 
conserving pseudopotentials [|ri|] in the Kleinman-Bylander form [|T^, treating the d-wave 
part as the local component. The exchange- correlation potential is computed using the local- 
density approximation (LDA), using the Ceper ley- Alder result for the exchange-correlation 
energy, as parametrized by Perdew and Zunger |]T3. The details of the code can be found 
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in the literature njTS 



Our liquid-state molecular dynamics (MD) simulations were carried out at a temperature 
T = 1600 K just above the melting point of GaAs (T = 1515 K). We have considered 
five concentrations of Ga^^Asi-a;: x = 0.2, 0.4, 0.5, 0.6, and 0.8. In each case, we used a 
cubic 64-atom supercell with periodic boundary conditions. For this size cell, the actual 
numbers of Ga atoms in the five samples were 13, 26, 32, 38, and 51. (Since the minority 
component in the 20% and 80% samples, the statistics will be quite poor for them and the 
results should be given greatest significance for the majority components in those cases.) 
The atomic densities for the five concentrations were obtained from the measured density 



of £-GaAs at this temperature |jT^, together with Vegard's Law (i. e. linear interpolation 
of atomic volumes) |jl8[ . We use a 10-Ry cutoff for the energies of the plane waves included 
in the wave function expansion, and F-point sampling for the supercell Brillouin zone. For 
the electronic structure, we used Fermi-surface broadening corresponding to an electronic 
subsystem temperature of /c^T^' = O.leV. In calculating the electronic wave functions, at 
each concentration we include eight empty bands. We control the ionic temperature using 
the Nose-Hoover thermostat [^]. The equations of motion are integrated by means of the 
Verlet algorithm, using an ionic time step of 125 a.u. (~ 3 fs). For each concentration, the 
samples were equilibrated for about 0.2ps, following which simulations were carried out for 
more than 3ps. 

These simulations start from an initial configuration which, for the stoichiometric liquid, 
was generated from a classical molecular-dynamics (CMD) simulation using potentials of the 
Stillinger- Weber form. Although our particular potentials were originally derived for liquid 
Ge^^Sii-a; alloys, they should give a reasonable starting configuration for the Gai-^As^; alloys 
which we study. We use these potentials but with Ga and As masses rather than those of Si 
and Ge. This model system was melted using classical molecular dynamics at high temper- 
atures, then cooled down gradually to the temperature of interest. The configuration thus 
obtained from CMD was equilibrated for 0.2 ps (about double the expected relaxation time 
for this system). For the nonstoichiometric liquids, we used the same starting configuration 



as in the stoichiometric case, but with Ga atoms randomly substituted for the appropriate 
number of As atoms (or vice versa), so as to give the correct concentration. 

III. RESULTS 

Fig. 1 shows the three partial pair correlation functions gcaCaij), gAsAsij) and gcaAs for 
liquid Ga^Asi-x. at the concentration x = 0.5 and a temperature T = 1600 K. A number 
of features deserve mention. First, at x = 0.5, the principal peaks of all three partial pair 
correlation functions occur at about the same separation, namely 2.5A. This indicates the 
non-ionic character of the bonds; in ionic liquids the partial pair correlation functions for like 
atoms are out of phase with the corresponding function for unlike atoms [^. Our results 
for ga/3(r) are in close agreement with the recent calculation by Godlevsky a/ an earlier 
calculation by Zhang et al ^ also gives the same features although these authors get a much 
stronger principal peak for gGayis(r) than is found in either our calculations of that of Ref. 

While the principal peaks are in phase, there are some differences among the partial 
g(r)'s. For example, at x = 0.5, gAsAs(r) has a slightly higher and narrower first peak, and 
a stronger second peak, than gGaGa(r), while gcaCaij) has a broad first peak and no obvious 
peak beyond that. We have also calculated the coordination numbers for the first shell of 
neighbors, defined as the integral of 47rpr^g(r) from zero out to the minimum after the first 
maximum in g(r). [Here g(r) is the tote/ pair correlation function, which does not distinguish 
between the two species, normalized so that it approaches unity at large r.] The coordination 
we calculate in this way at T = 1600 K is 5.8, in good agreement with the experimental 
estimate of 5.5 ± 0.5. Note also that this value is larger than the diamond-structure value of 
4 but significantly smaller than the value expected in a close-packed liquid, which would be 
in the range of 9 or 10. This value indicates the persistence of covalent bonding in ^-GaAs. 

Fig. 2 shows gGaGa(j) sX X = 0.8 and gAsAsij) at a; = 0.2. The latter shows more short- 
range order than the former - specifically, a sharper main peak and a broad second peak, 
rather than a single broad principal peak. We believe that there are several causes for these 
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differences. First, pure Ga has a much lower melting temperature than either Gao.sAso.s or 
As. Thus, at the same temperature of 1600 K, we expect less short-range order for gGaGa(r) 
at X = 0.8 than for g^s^s(r) at x = 0.2, as seen in our calculations. In addition, we expect 
the X = 0.2 sample to show some residue of the complex local structure seen in pure i-As, 
which is, in turn, quite similar to that of the crystalline phase P,p^ . Indeed, our calculated 
SAsAsi'f') at x = 0.2 has some of the same features as those in the calculated g(r) for pure 
i-As [|]. However, there are also some observable differences, which may be related to the 



fact that, experimentally, pure i-As is semiconducting p2[ while £-Gao,2Aso.8 is calculated to 
be metallic (as shown later when we calculate the electronic structure). First, if we integrate 
gAsAs(r) for the x = 0.2 sample out to the first minimum beyond the main peak, we obtain 
a coordination number of 3.2 As atoms for the first shell of nearest neighbors surrounding 
an As atom. This is slightly larger than the value 3 reported for pure £-As (at T=1150 K) 
experimentally and in previous calculations |]8|j23[|. A more important difference is that pure 
^-As has a sharp second peak in g(r) [^, whereas for gAsAs(r) in our simulations this second 
peak is very broad. While some of this broadening may result from the higher temperature 
(T = IGOOK) in our simulations for x = 0.2, these results suggest that the local structure 
seen in i-As is preserved only out to the first shell of neighbors in £-Gao.2Aso.8- 

On the other hand, the features we see for gG'aGa(i') at x = 0.8 are qualitatively similar and 
consistent with those seen for liquid Ga at lower temperatures Using the procedure 

indicated above, we get a Ga coordination number of 5.8 for the first shell. If we integrate the 
total g(r) up to the first shell (i. e., including both As and Ga neighbors of Ga), we obtain 
a Ga coordination number of 6.9. While these values are smaller than the coordination 



number of 9.0 reported for pure i-Ga p6[, we attribute the difference to the lowering of the 
first peak due to the higher temperature of our simulations (T = 1600 K) compared to those 
for pure i-Ga {T = 702 K and T = 982 K) §. 

Figs. 3-5 show information about the various partial alloy structure factors at the same 
temperature. The partial structure factors Sij(k) are defined in one of the standard ways 

MM- 
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where i and j denote the two components of the binary alloy. Fig. 3 shows the calculated 
total structure factor S(k), as weighted by neutron scattering factors at a concentration 
X = 0.5. S(k) is defined by 

_ k'Suik) + 2kb,S,,ik) + b,'S,,ik) 

where hi and bj are the corresponding experimental neutron-scattering lengths ( bca = 7.2 
and bAs = 6.7) For comparison, we also show in Fig. 3 the quantity S(k) as measured by 



neutron diffraction As can be seen from the Figure, the two agree quite well. In 

particular, the calculations convincingly reproduce the experimentally observed shoulder on 
the high- A; side of the principal peak in S(k). 

Fig. 4 shows the three partial alloy structure factors ScaCaO^), SasAsO^), and SGaAs(k) 
at X = 0.5 and T = 1600K. The structure factor for like pairs is always positive, with a 
conspicuous first peak, while the structure factor between opposite pairs is negative for small 
k, becoming positive at k values corresponding to the peaks in the other two partial structure 
factors. It is of interest to compare these results to those found in other model calculations. 
For a mixture of hard spheres of packing fraction 0.45 (characteristic of the liquid near 
melting) and ratio of hard sphere diameters of 0.9 [^], the cross-correlation function is 
negative at small k and has a peak near that of the two same-species functions, as in our 
calculations. By contrast, for liquid NaCl, which is strongly ionic), the cross-correlation 
function has a strong negative peak at the same k as the peaks of the same-species partial 
structure factors. Thus, our results are more similar to the hard-sphere structure factors, 
suggesting that liquid GaAs is at most only weakly ionic. 

Finally, Fig. 5 shows ScaGaiX) at x = 0.8 and S^s^s(k) at x = 0.2. Once again, like 
the real-space correlation functions at the same concentrations, the Ga-Ga structure factors 
show slightly less correlation (i. e., a lower principal peak and a less conspicuous second peak) 
than do the corresponding As-As structure factors. One possible reason for this behavior, as 
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for the corresponding real-space correlation functions, is that the x = 0.8 sample is further 
from melting than is the x = 0.2 liquid. 

SA<iAs(k) for X = 0.2 shows characteristic differences from that for pure i-As, which are 



analogous to those discussed earlier for the g(r)'s. Specifically, S(k) for pure i-As |^ has a 
split principal peak with maxima at A; = 2.45 and k = 3.74 A~^. By contrast, SAsAs(k) 
for X = 0.2 has a peak at 2.5 A, but the second peak is reduced to only a shoulder at about 
3.5 A. The fact that the second of the split peaks is smoothed to a shoulder at x = 0.2 
indicates a change in the local structure which is also reflected in the reduction of the second 
peak of partial g(r) as noted earlier. 

We next discuss the x = 0.8 sample, comparing our results with those for pure i-Ga at 
lower temperatures P,|26[|. The partial structure factor show the same qualitative features 



as the structure factors for the pure liquid, but there are some quantitative differences. Our 
results for x = 0.8 (at 1600 K) show a first peak in ScaCaiX) with a maximum of only about 
1.3, which is lower than the experimental one seen in pure Ga at 959 K (~ 1.7). We attribute 
the lowering of the first peak to the increased temperature, as seen experimentally in most 
liquid metals. 

Further information about the short-range order in the liquid alloy may be obtained from 
bond angle distribution functions, shown in Figs. 6 and 7. These functions are defined by 
analogy with our previous work in liquid Ge p|JTO[|. Namely, one considers a group of three 
atoms. Of these, one is denoted as the central atom; the other two atoms (denoted as "side 
atoms"), with the central atom, define a bond angle 9. g3{0,rc) is the distribution of bond 
angles formed by all such groups of three atoms, such that both the side atoms lie within 
a cutoff distance Tc of the central atom. Fig. 6 shows g^{9,rc) for the Ga-Ga-Ga angles at 
X = 0.8, and for the As- As- As angles at a; = 0.2, each for two different choices of the cutoff 
radius r^. Fig. 7 shows the same functions for Ga-Ga-Ga and for As-As-As at x = 0.5. 

The peaks in these distribution functions give hints about about the short-range bond- 
order in the liquids. For example, a peak near 9 = 60° corresponds to a relatively close- 
packed arrangement of the corresponding atomic group, with many nearest neighbors. In 
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contrast, a peak near 100" indicates a more tetrahedral structure, typical of covalent bonding. 
Thus, the upper part of Figs. 6 and 7 suggest that the Ga ions form a rather close-packed 
arrangement at x = 0.8 and x = 0.5, since there is a strong peak near 60" for both values 
of the cutoff Tc- By contrast, the lower parts of Figs. 6 and 7 suggest that, as expected, the 
As atoms have a more open arrangement at x = 0.2, since there are two peaks in at this 
concentration: a side peak near 50°, and another noticeable peak near 97°. We also observe 
the peak at 97° for As atoms at x = 0.5, however it is broader and less pronounced than at 
X = 0.2. In both Figs. 6 and 7, the Ga-Ga-Ga bond angle distribution depends little on Tc. 

At X = 0.2, the As- As- As distributions show noticeable 50 — 60° peaks only at the larger 
cutoff values. At the smaller cutoff radius, the 50 — 60° peak is missing. Thus, at short 
distances, the As clusters tend to maintain the local version of the structure they have in 
the pure liquid phase (and the crystalline phase) which shows a strong peak at 97°, but at 
larger cutoffs the local structure differs from that of pure £-As. We have made the same 
observation in connection with gAsAs(r) at the same concentrations, and with the shoulder 
in SAsAs(k). The 97° peak implies some tetrahedral order persisting to x = 0.2, though this 
peak is less pronounced than in pure i-As 0. 

We have also calculated the electronic properties of i-Ga^Asi^x- We calculate the single- 
particle electronic density of states N{E) in the standard way, by using the expression 

N{E)= J2w^g{E-E^). (3) 

In this expression i^k denotes one of the energy eigenvalues of the single-particle wave func- 
tions at a particular k point within the supercell Brillouin zone, Wk is the weight of that 
k point (defined below), and g(E) is a Gaussian smoothing function of width a = 0.2eV. 
Our calculation is carried out by sampling the supercell Brillouin zone at eight special k 
points, using the same choice of special points and weights as that of Holender et al in their 
calculations for pure £-Ga 0, which have been well tested and found to be an adequate rep- 
resentation of the supercell Brillouin zone. This choice is, however, convenient rather than 
unique; we expect that other choice having the same number of k points would have given 
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similar results, as has been found in pure £-Ga. For each k-point we include 40 conduction 
band states, and for each concentration, we obtain our final results by averaging over twelve 
representative liquid state configurations. 

The resulting calculated density of states N(E) is shown for the four concentrations 
X = 0.2,0.4,0.6, and 0.8 in Fig. 8. [We do not show our calculated N(E) for x = 0.5, but 
it interpolates smoothly between x = 0.4 and x = O.6.] As in our previous studies, the 
alloy has a clearly metallic density of states for all concentrations x. However, just as in our 
previous results for Gaa;Gei_2^, the density of states becomes more and more free-electron 
like as the concentration x of the metallic component (Ga in this case) increases. Pure £-As 
is semiconducting and has been calculated to have deep minima in the electronic density 
of states at the Fermi energy (Ep), and also at an energy of ~ -7 eV (measured from Ep) 

which separates the s and p bands. Liquid Ga, on the other hand, has an almost free- 
electron like density of states. We see these features reflected in our simulations; as the 
figures indicate, the electronic density of states (DOS) has a pseudogap in the As rich phase 
which progressively fills up as the Ga concentration is increased, so that for x = 0.8 it is 
hardly noticeable. But even at low Ga concentration {x = 0.2), there is no minimum in 
the density of states at the Fermi energy. As for the pseudogap, we find that its position 
changes monotonically to lower energies (relative to Ep) with increasing Ga concentration. 
A similar pseudogap is reported in calculations for pure £-Ge for the same reason, i. 

e., a partial separation between s-like and p-like bands. 

We have also computed the frequency- dependent electrical conductivity (y{u>) for x = 0.2, 
0.4, 0.5, 0.6, and 0.8, at frequencies ranging up to 2eV. cr^uj) is given by the standard Kubo- 
Greenwood expression pO| 
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where t/'j and Tjjj are the single particle Kohn-Sham wave functions with Fermi occupancies 
fi and fj and energy eigenvalues Ei and Ej. Once again, we calculate the conductivity 
using the same set of eight special k points used for N(E), and averaging over the same 
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twelve representative ionic configurations, including 40 conducting band states for each k. 
While there is some inaccuracy in using differences in Kohn-Sham eigenvalues, as we do 
in this equation, it is difficult to do better at present in liquid metals and semiconductors, 
and this same approach does give good agreement with experiment in a number of other 
liquid metals, especially at zero frequency but also for stoichiometric ^-GaAs ^ at finite 
frequencies. 

The frequency-dependent conductivity is shown in Fig. 9 for x = 0.2 and x = 0.8, and 
its calculated zero-frequency limit is listed in Table I for the five concentrations x = 0.2, 0.4, 
0.5, 0.6 and 0.8, all at T = IGOOiT. 

Several features of the conductivity graphs, and of the tabulated d. c. limits, deserve 
mention. First, the calculated value of the d. c. conductivity at x = 0.5 is very close to the 
experimental value: 0.84 x 10^ ohm^^cm^^, compared to the experimental value of 0.79 x 10^ 
ohm~^cm^^ Secondly, the conductivity has a weak minimum near x = 0.5. This is 

consistent with expectations based on second-order perturbation theory which would 
predict that alloy scattering (due to concentration fluctuations) would be a maximum near 
X = 0.5 Thirdly, the frequency- dependence of cr{uj) becomes more metallic as x increases. 
Specifically, at x = 0.8, (j{uj) clearly decreases with increasing frequency, characteristic of a 
Drude metal, while for the highest As concentration (x = 0.2), the conductivity is nearly 
frequency-independent in the range of calculation. This behavior is closer to the nonmetallic 
behavior in which the conductivity increases with increasing frequency. 

Finally, we have computed one important atomic transport coefficient, namely, the 
atomic self-diffusion coefficients Da for the majority species in the two liquids Gao.2Aso.8 
and Gao.8Aso.2- In both cases, the Da's can be extracted from a plot of the mean-square 
atomic displacement versus time, which approaches a straight line in the limit of large time. 
The expression is 

_ (|H.W-H.(0)P>^ (5) 

where Ri(t) denotes the position of an ion of species i at time t, and the triangular brackets 
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denote an average over all atoms of species i and over initial times. Plots of the mean-square 
atomic displacement as a function of time are shown in Fig. 10 for both types of atoms at 
X = 0.2 and 0.8 as indicated above; iFiom these we obtain Dca = 2.1 x 10~^ cm^/sec at 
X = 0.8; Das = 1.7 X 10"^ cmVsec at x = 0.2. 

The value for D^s at x = 0.2 is about three times larger than that obtained by Li [^] for 
pure £-As [D ~ 0.6 x 10~^ cm^/sec). This difference is probably due to a combination of 
several factors. First, the calculations by Li are carried out at a temperature ~ 450 K lower 
than ours. Secondly, pure £-As seems to have more covalent bonding than Gao,2Aso.8, which 
probably impedes atomic motion, giving a lower atomic diffusion coefficient for As. Thus, 
in short, this behavior seems to be consistent with the rest of our picture, which is that 
^-Ga^^Asi-j; rapidly acquires metallic conductivity, and corresponding structural properties 
for X as small as 0.2. 

IV. DISCUSSION AND CONCLUSIONS 

The most striking results of these calculations is that Ga^^Asi-a; remains metallic at all 
concentrations between x = 0.2 and x = 0.8, including the As rich value x = 0.2. Thus 
even low Ga concentrations are sufficient to render i-Ga^Asi^^ metallic (recall that pure 
^-As is semiconducting). This is reflected in the electronic density of states which shows no 
minimum at the Fermi energy for all concentrations studied. The liquid structure is also 
consistent with metallic behavior at all concentrations between x = 0.2 and 0.8, although 
there are some clear deviations from the behavior seen in simple metallic alloys. Specifically, 
although the coordination number at all concentrations is larger than the value of four that 
might be expected in a predominantly covalent liquid, it is still smaller than that of a typical 
close-packed hard-sphere mixture. At a: = 0.2 and 0.8, the pair correlation functions and 
structure factors resemble those of the corresponding pure liquids, except that the split first 
peak in S{k) of i-As becomes a single peak with a weak shoulder in As. 

It has been suggested that semiconducting properties persist in a liquid only if the 
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liquid has the same short-range order as the crystalhne phase. In the alloys we study, the 
liquid state has a short range order which is distinctly different from the solid. We also find 
that all these alloys are metallic, in agreement with the suggestion of Ref. |3^. By contrast, 
again in agreement with the picture advanced in another compound semiconductor, 
stoichiometric CdTe, appears to preserve the crystalline short-range order, and also to retain 
its semiconducting characteristics in the liquid state. Such semiconducting behavior was 
indeed found in the ab initio calculation by Godlevsky et a/ |Q. It would be of interest to 
extend their calculations off stoichiometry, where metallic behavior is likely. 

The other characteristics of ^-Gaj^Asi^i; reflect its fundamentally liquid-metallic charac- 
ter. For example, the resistivity is predicted to exhibit typical Faber-Ziman behavior: a 
weak positive deviation from a linear concentration dependence near stoichiometry, which is 
caused by alloy scattering. The calculated values of the partial atomic diffusion coefficients 
are comparable to the diffusion coefficients of those other liquid semiconductors which are 
metallic in their liquid states. We find no evidence of a strong reduction in this value be- 
cause of formation of clusters near stoichiometry; such cluster formation is not expected for 
Ga^Asi„2: because of the small electronegativity differences between the two species. One 
might speculate that, in other liquid semiconductors, such as stoichiometric i-CdTe, which 
remain poorly conducting in the liquid state, the local structure is much more ionic near 
melting and the atomic diffusion coefficients are correspondingly lower. 

In summary, our calculations show that Gaa;Asi_a; is a reasonable metal at all concentra- 
tions between x = 0.2 and x = 0.8. In particular, there is no evidence of strong compound 
formation in the liquid state near x = 0.5. The electrical conductivity shows a concentration 
dependence typical of a liquid metallic alloy, with evidence of weak scattering from concen- 
tration fluctuations which reaches a maximum near x = 0.5. The electronic density of states 
shows no minimum at the Fermi energy; instead, it has a pseudogap between s-like and 
p-like occupied state which persists at all concentrations, though it is considerably weaker 
in the Ga-rich alloys. The atomic diffusion coefficient is calculated to be similar to that 
of other liquid semiconductors which are metallic in their liquid state. Finally, the liquid 
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structure shows some indications of deviation from the behavior of simple hquid metal alloys 



25| . The principal evidence of deviation from behavior characteristic of a simple liquid alloy 



formation comes from the calculated pair correlation functions, structure factors, and bond 
angle distribution functions, all of which show some weak indications of departure from 
close-packed behavior: smaller coordination functions than in typical hard-sphere liquids, 
and a weak residue of tetrahedral bonding. 
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TABLES 

TABLE I. D. c. conductivity at the five concentrations, obtained by extrapolating low frequency 
a. c. conductivity results. 



Concentration (%) 


20 


40 


50 


60 


80 


Gdc (10^ ohm"^ cm~^) 


0.93 


0.9 


0.84 


0.91 


1.1 
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FIGURE CAPTIONS 



1. (a) Pair correlation functions gGaGa(j), gAsAsij), and gcaAsij) for £-Ga^,Asi_^. at a 
temperature of 1600 K and concentration x = 0.5. The graphs are vertically offset by 
one unit each for clarity. 

2. (top) gGaGa(r) for X = 0.8; (bottom) gAsAsir) at x = 0.2 and T = 1600 K. 

3. Full line: Calculated neutron structure factor S(k) = X)i^j=Ga%(k)Sij(k), where fij{k) 
is the neutron form factor, for x = 0.5. Open circles: measured S(k) as obtained by 
neutron diffraction [ref. p7| ]. 

4. Partial structure factors ScaGaO^), ^AsAsiX), and SGayis(k) for x = 0.5. S^syis(k) is 
vertically offset for clarity. 

5. Calculated partial structure factors (top) SGaGa(k) at x = 0.8; (bottom) SAsyis(k) at x 
= 0.2. 

6. Calculated bond angle distribution functions g{9, Vc) for (top) groups of three Ga atoms 
at a concentration x = 0.8, and (bottom) groups of three As atoms at a concentration 
X = 0.2, both for two different cutoff radii, Tc = 3AA and Tc = 3.8A, as defined in the 
text. 

7. Same as for Fig. 6 but for x = 0.5 and (top) gGaGaGa(^' '"c); (bottom) gAsAsAs^^' ^c)j 

8. Single-particle electronic density of states N{E) for Ga^^Asi-^: at x = 0.2, 0.4, 0.6, and 
0.8. The Fermi energy in each case is shown as a dashed vertical line. 

9. Calculated electrical conductivity cr(u;) for Ga^Asi-^: at a; = 0.2 and x = 0.8. 

10. Calculated mean-square displacement (|Rj(0) — Ri(it)p), plotted as a function of t at 
X = 0.2 for As and at x = 0.8 for Ga. 
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FIG. 4. 
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FIG. 5. 
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FIG. 6. 
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FIG. 7. 
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FIG. 9. 
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